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Abstract 

The Great Lakes of East Africa are among the world's most important freshwater ecosystems. Despite their importance in 
providing vital resources and ecosystem services, the impact of regional and global environmental drivers on this lacustrine 
system remains only partially understood. We make a systematic comparison of the dynamics of the bio-optical and thermal 
properties of thirteen of the largest African lakes between 2002 and 201 1. Lake surface temperatures had a positive trend in 
all Great Lakes outside the latitude of 0° to 8° south, while the dynamics of those lakes within this latitude range were highly 
sensitive to global inter-annual climate drivers (i.e. El Nino Southern Oscillation). Lake surface temperature dynamics in 
nearly all lakes were found to be sensitive to the latitudinal position of the Inter Tropical Convergence Zone. Phytoplankton 
dynamics varied considerably between lakes, with increasing and decreasing trends. Intra-lake differences in both surface 
temperature and phytoplankton dynamics occurred for many of the larger lakes. This inter-comparison of bio-optical and 
thermal dynamics provides new insights into the response of these ecosystems to global and regional drivers. 
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Introduction 

The Great Lakes of East Africa are among the world's most 
important aquatic ecosystems from the point of view of ecosystem 
services, biodiversity and carbon cycling [1]. They are highly 
heterogeneous, presenting an array of hydrological and biological 
characteristics, ranging from shallow eutrophic water bodies (e.g. 
Lake Victoria) to deep oligotrophic lakes (e.g. Lake Tanganyika). 
Recent studies have demonstrated linkages between climate 
variability and productivity of several African lakes [2,3,4]. This 
is partially due to the importance of direct rainfall and evaporation 
in the water balances in the larger lakes (e.g. Lakes Tanganyika, 
Malawi and Victoria) and river flow in the smaller lakes (e.g. Lake 
Turkana) as well as climate driven lake stratification cycles [5,6,7]. 
Inter-annual variations in regional and global climate patterns 
have a strong influence on monsoon dynamics, in particular the 
Inter Tropical Convergence Zone (ITCZ) [8] . Local and regional 
impacts also play an important role on the chemical and biological 
dynamics of many of these lakes [9,10,1 1]. Remote measurements 
have been used to explore the spatio-temporal variability of the 



temperature and bio-optical properties in several of these lakes 
[12,13,14]. However, a systematic inter-lake comparison using a 
common methodology has yet to be performed. There are 
multiple reasons for this, geographical (e.g. altitude) and seasonal 
differences in atmospheric optical conditions and seasonal 
differences in aquatic optical properties. Most importantly, there 
are limited data for lake specific algorithm development and 
validation while standard bio-optical algorithms have been 
predominately developed for phytoplankton dominated water 
bodies found at sea level [15]. In the present study, we make a 
comparative analysis of the thermal and bio-optical dynamics of 
thirteen lakes lying in or between the two branches of the Great 
Rift in Eastern-Central Africa; Lakes Albert, Chilwa, Edward, 
Kivu, Kyoga, Malawi, Mweru, Naivasha, Rukwa, Tana, Tanga- 
nyika, Turkana and Victoria (Figure 1). Using site-specific seasonal 
smoothing of the time series, we examine their decadal trends. 
These trends are also compared to interannual climate drivers, 
allowing us to explore commonalities between lakes and potential 
drivers of change. 
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Figure 1. Lakes and lake sections included in the present 
analysis of the African Great Lakes. 

doi:1 0.1 371 /journal.pone.0093656.g001 

Materials and Methods 

Monthly estimates of chlorophyll-a concentrations and lake 
surface temperatures (night) were obtained as Level 3 data from 
the MODIS AQUA mission dataset from the Goddard Earth 
Sciences Data and Information Services Centre (http://disc.sci. 
gsfc.nasa.gov/giovanni, Ocean Color Radiometry Online Visual- 
ization and Analysis Monthly Data), following Reprocessing 1.1 
[16] and covering the period from July 2002 to July 2011. 

Lake sections (64 km 2 -128 km 2 ) were delineated using MODIS 
AQUA data at 4 km spatial resolution for the bio-optical estimates 
(MAMO_CHLO_4 km). Surface temperature estimates of the 
same lake section were obtained using MODIS AQUA data at 
9 km spatial resolution (MAMO_NSST_9 km). Larger lakes were 
divided into separate lake sections based on the large scale patterns 
identified in previous studies or where lake dimension suggested 
such patterns could be possible [17,14]. Possible influences of 
coastal waters and edge effects were minimized in the large lakes 
by selecting lake sections at least 10 km from the lake border. 

To avoid errors due to geographical differences between the 
optical conditions (atmospheric and aquatic) of the African lakes 
and those of standard MODIS processing, each record category 
was converted to a vector of anomalies by subtracting the lake 
wide average and dividing by the lake wide standard deviation. In 
such a manner, consistent spatial differences in optical conditions 
between lake sections would not influence the analysis of lake 
dynamics [17]. 

To remove the influence of seasonal variability in optical 
conditions (atmospheric and aquatic) between lakes and lake 
sections, each vector underwent a seasonal smoothing by using a 
yearly centred moving average with a length of 12 (months). The 



resulting trend vector for each parameter was then subtracted 
from the original data and the differences were averaged across all 
contemporary (monthly) observations to produce a seasonal index. 
The trend vector and seasonal indices were subtracted from the 
data to give a residual component which was checked for normal 
distribution. 

The resulting trend vectors (SST for the trend vector of lake 
surface temperature anomalies, CHLa for the trend vector of 
chlorophyll a anomalies) were independent of seasonal and intra- 
lake spatial variations in optical conditions (atmospheric and 
aquatic). These vectors of seasonally smoothed data anomalies 
were then used to compare dynamics between lakes and lake 
sections to determine commonalities in inter-annual trends. 

Seasonally smoothed trend vectors of meridional wind speed 
anomalies (WIND) and precipitation rate anomalies (RAIN) for 
each lake section were constructed using NCEP/NCAR data 
[18,19] (http:/ /www.esrl.noaa.gov/ psd/data/gridded/ data.ncep. 
reanalysis.html). Likewise, a trend vector of the Multivariate El- 
Nino Southern Oscillation Index (MEI) was also constructed [20] 
(http://www.esrl.noaa.gov/psd/enso/mei/). To examine the role 
of the ITCZ, we used historical dekadal (10-day) data for the 
African position of the ITCZ from April 2003 to July 2011 
provided by the NOAA CPC Africa Desk/International Desk's 
Team (www.cpc.ncep.noaa.gov/products/fews/ITCZ/itcz) [21]. 
Monthly averages of dekadal positional data for longitudes of 30 
and 35 were averaged and compared to SST and CHLa 
anomalies. 

Factor analysis was utilised to identify common temporal co- 
variation between lake sections. This analysis was based on a linear 
combination of variables which contained the maximum variance. 
The reduced numbers of factors (eigenfunctions) is orthogonal 
(uncorrelated) while the eigenvalue of each factor was used to 
determine the fraction of the variability explained. In this manner, 
the temporal dynamics of the 27 lake sections was reduced to a 
smaller number of dominant factors. Common temporal modes 
were identified by determining the correlation between each lake 
section and the dominant factors for each variable. Lake sections 
with similar temporal dynamics were correlated to the same factor. 
It should be noted that this approach assumes that the temporal 
variability of each variable dataset is due to common forcing 
agents. 

Decadal trends and correlations between drivers and anomaly 
vectors were determined by using linear model with a monthly 
time scale where direct correlations (increasing) had positive 
Pearson correlation coefficients (r) and indirect correlations 
(decreasing) had negative coefficients. 

Results 

The SST trend vectors of each lake section (27), analyzed by 
factor analysis, had five significant temporal modes which 
accounted for 92.5% of the total variability of the dataset 
(60.2%, 14.0%, 9.0%, 5.6%, 3.8% respectively). Sections of the 
same lake most often followed the same temporal mode (Table 1). 

The SST trend vector for each lake section from 2002-2011, 
examined using a linear model and decadal trends, was identified 
as either positive (increasing), negative (decreasing) or not 
significant at p>0.01 (Table 1, Figure SI in File SI). The decadal 
trend for the SST vector was positive for Lakes Albert, Chilwa, 
Edward, Kyoga, Malawi (all sections), Mweru (both sections), 
Rukwa, Tanganyika (south and southcentre), and Turkana (both 
sections), indicating an overall warming of the lake surface waters. 
A negative trend was observed for Lake Tana. The highest inter- 
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annual variability (lowest r) was found in those lakes south of the 
equator (Lakes Victoria, Kivu, Tanganyika (north)). 

Using the CHLa trend vectors, factor analysis identified five 
significant temporal modes that accounted for 80.8% of the 
variability of the dataset (27.0%, 24.8%, 13.5%, 8.6%, and 6.9%, 
respectively). Differences in the inter-annual trend between 
sections of the same lake were evident in Lakes Tanganyika and 
Victoria. Positive trends (p<0.01) in CHLa occurred in Lakes 
Albert, Edward, Malawi (all sections), Mweru (both sections) and 
Victoria (centre) (Table 1, Figure SI in File SI). Negative trends 
were evident in Lakes Kyoga, Tanganyika (north and north 
centre) and Turkana. The highest inter-annual variability was in 
the smaller lakes as well as the south section of Lakes Tanganyika 
and Malawi. 

Factor analysis of WIND and RAIN identified three and four 
dominant modes that accounted for 89.9% and 95.2% of the 
variability of each dataset. Negative trends in WIND were 
observed in most lakes at or above the equator, while positive 
trends occurred in lakes below the equator. Positive trends in 
RAIN were observed in all lakes except Lakes Malawi and Mweru 
(negative) and Lake Tanganyika (not significant). 

The multivariate ENSO index (MEI), used to compare inter- 
annual global climate variations with the SST and CHLa 
dynamics in each lake section (Table 1, Figure S2 in File SI), 
showed a significant positive correlation between MEI and SST 
(p<0.01) in sections of Lakes Victoria, Tanganyika, Kivu and 
Mweru, while a negative correlation for Lake Chilwa. Correlations 
between CHLa trends and MEI were both positive and negative 
(Table 1) with intra-lake differences. The decadal trend of MEI 
during the study period was negative and included the 2007-2009 
and 2010-2011 La Nina events and the 2009-2010 El Nino 
events. WIND was positively correlated to MEI in most lakes while 
RAIN was negatively correlated. 

The comparison of SST anomalies to ITCZ positional data 
(Table 1) showed a significant negative correlation between ITCZ 
and SST (p<0.01) in all lake sections with the exception of Lakes 
Tana and Turkana (north). The correlation increased at lakes with 
latitudes below the equator. Correlations between CHLa anom- 
alies and ITCZ positional data were positive in several lakes, most 
significantly in Lakes Victoria and Mweru (Figure S2 in File SI). 

Discussion 

The decadal warming observed in both large and small lakes 
results from inter-annual variability in radiative balances and 
latent heat exchange, as well as indirectly from changes in lake 
mixing and lake water balances. Interestingly, significant warming 
trends were limited to those lakes which displayed the lowest 
correlation with MEI. SST trends in Lakes Victoria, Tanganyika 
(north, north centre and centre) and Kivu (both sections) showed a 
clear sensitivity to MEI, which led to higher inter-annual 
fluctuations but no significant decadal trend. The spatial pattern 
of the MEI-SST relationship confirms that the lakes between 0° to 
8° south are most sensitive to El-Nino effect on the ITCZ and the 
related weakening of the monsoon wind flow over East Africa [23]. 
Our observations suggest that the effect of El-Nino phenomenon 
on surface heat balances can be recurrent and significant for many 
lakes in Tropical East Africa [14]. The high correlation between 
the SST dynamics and the latitudinal movement of the ITCZ 
(Table 1) confirms the influence of climate drivers on the dynamics 
of the African lakes [21,22]. 

Trends in CHLa are influenced by inter-annual variability in 
nutrient availability and optical conditions related to light 
limitation. Both factors are related to long term trends in local 



and regional conditions, as well as climate. The Lake Albert CHLa 
trend had a high negative correlation with MEI while Lake 
Turkana CHLa (north and south) was positively correlated. 
Opposing correlations between CHLa trends and MEI in the 
north and south of Lake Tanganyika were found. The influence of 
the latitudinal movement of the ITCZ was evident in several lakes 
(Lakes Kivu, Mweru and Victoria) and individual sections of Lakes 
Malawi and Tanganyika. 

Decadal trends in CHLa in Lake Tanganyika also showed 
latitudinal differences, while all sections of Lake Malawi demon- 
strated a positive trend in CHLa over the study period, indicating 
a general increase in phytoplankton biomass. The centre section of 
Lake Victoria showed a significant positive trend, while the south 
section showed a similarly positive but less significant (p = 0.025) 
decadal trend. 

CHLa dynamics are influenced by lake hydrodynamics through 
variations in seasonal cycles in mixing depth and nutrient 
upwelling. While SST dynamics are based on thermal emission 
from surface waters, surface temperature anomalies have been 
used to track changes in lake hydrodynamics [12,14]. Lake 
Tanganyika CHLa (south, southcentre and centre) was negatively 
correlated to SST. Such a negative relationship would indicate 
that upwelling nutrient rich waters have a positive impact on 
phytoplankton productivity [12]. A similar negative relationship 
was found for the south section of Lakes Kivu. Lake Victoria 
behaved in an opposite manner, as CHLa and SST were 
negatively correlated in the north section only, with a weak, but 
positive relationship in the south. A positive relationship between 
CHLa and SST could indicate a positive effect of lake stratification 
on phytoplankton growth (e.g. due to light limited growth). The 
CHLa trend in the north and centre sections of Lake Malawi 
showed a positive correlation to SST, as did correlations between 
SST and CHLa in Lakes Edward, Mweru (north) and Turkana 
(north). 

Intra-lake Variability in Trend Vectors 

In the present dataset, intra-lake differences in trend vectors of 
SST (Table 1) were significant in Lakes Malawi, Tanganyika, 
Victoria and Albert (p<0.001) and to a lesser degree in Lakes 
Mweru (p = 0.01). In Lakes Victoria, Malawi and Mweru, the 
northern SST trend vectors were continuously higher than the 
southern lake sections, indicating a stable latitudinal SST 
difference in each lake. Lake Albert SST trend vectors were 
highest in the south section. It should be noted that seasonal 
latitudinal switches in SST may also be present but could remain 
undetected using anomaly trend vectors. Large scale latitudinal 
differences in surface temperature have been associated to 
differences in the exposure to the monsoon winds and the 
duration of the windy season (e.g. Lake Victoria [14]). On the 
other hand, the highest SST trend vectors in Lake Tanganyika 
were found in the north and south sections while the lowest SST 
vector occurred in the centre section. Large scale latitudinal 
circulation has been reported in Lake Tanganyika [24,25]. 

Intra-lake differences in CHLa trends were significant in most 
lakes. In Lakes Tanganyika and Malawi, the CHLa trend vectors 
in the south were significantly higher than those of the northern 
sections. In Lake Tanganyika, latitudinal differences in chlorophyll 
a concentration have been confirmed in several studies 
[26,27,28,29]. Similarly, the south sections of Lake Malawi had 
been reported to have the highest productivity for most of the year 
[30,31]. In Lakes Albert and Mweru (Figure 2), the north section 
CHLa trend vectors were consistently higher than those of the 
south. Lake Kivu's CHLa dynamics showed less clear latitudinal 
differences, even though differences in phytoplankton composition 
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Figure 2. Trend vectors of A) CHLa for north and south sections of Lake Mweru, B) SST for north, centre and south sections of Lake 
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ENSO MEI. 

doi:1 0.1 371 /journal.pone.0093656.g002 



have been reported, with cyanobacteria more abundant in the 
north section while diatoms and cryptophytes being a higher 
proportion of total phytoplankton biomass in the south [32] . 

Conclusions 

The anomaly trend vectors of SST indicated that a decadal 
warming has occurred in all of the Great Lakes outside the latitude 
range of 0° to 8° south. Those lakes unaffected by the decadal 
warming trend were found to be more sensitive to inter-annual 
climate fluctuations related to the El Nino phenomena. It should 



be noted that the warming trends are limited to the study period, 
but the increased sensitivity of lakes in the latitudinal band of 0°— 
8° south to El Nino was clear. On the other hand, the sensitivity to 
the latitudinal position of the ITCZ was elevated in nearly all lakes 
south of the equator. 

The relative response of CHLa to changes in SST, MEI and 
ITCZ varied between lakes and between sections of the same lake. 
There was no clear regional trend towards eutrophication, as some 
lakes showed a negative decadal trend in CHLa, including sections 
of Lakes Tanganyika, Kyoga and Turkana. Lake Malawi presents 
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a positive decadal trend in CHLa and a positive relationship 
between SST and CHLa, indicating that this lake is moving to a 
higher trophic status, with a marked sensitivity to hydrodynamics. 

The regional and seasonal differences in the atmospheric and 
aquatic optical properties of the East African Great lakes, 
combined with limited field data have hampered a regional 
analysis. The present approach does not replace the need for 
extensive field data or for the development of lake specific 
algorithms to estimate aquatic and atmospheric parameters. 
However, the use of anomaly trend vectors allows for the 
identification of common lake trends and possible links between 
potential drivers and lake responses. 

Supporting Information 

File SI Figure SI. Decadal trends (2002-2011) in anomaly 
vectors of A) SST and B) CHLa in the African Great Lakes, blue 
circles for a positive trends, red circles for a negative trends, open 
circles for no significant trend (p>0.01). Figure S2. Correlations 
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